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In this paper we develop a new method that is different from Schwinger proper time method 
to deduce the fermion propagator with a constant external magnetic field. In the NJL model, 
we use this method to find out the gap equation at zero and non-zero temperature, and give the 
numerical results and phase diagram between magnetic field and temperature. Beside these, we also 
introduce current mass to study the susceptibilities, because there is a new parameter (the strength 
of external magnetic field) in this problem, corresponding this new parameter, we have defined a 
new susceptibility \b to compare with the other two susceptibilities y c (chiral susceptibility) and 
Xt (thermal susceptibility), and all of the three susceptibilities show than when current mass is not 
zero, the phase transition is a crossover, while for comparison, in the chiral limit, the susceptibilities 
show a second order phase transition. At last, we give out the critical coefficients of different 
susceptibilities in the chiral limit. 
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I. INTRODUCTION 

In recent years, it is widely believed that strong magnetic field could play an important role in astrophysics 0,1 
and high-energy physics [3H5[. Especially in quantum field theory, many papers have shown that external magnetic 
field can affect dynamical chiral symmetry breaking (DCSB), which is also known as ‘ Mag netic Catalysis’ j6HlO|. In 
the massless Nambu-Jona-Lasinio (NJL) model, with Schwinger proper time method [10 13] (one also can refer to 
Ref. [3] for another method), it is shown that in both 1 + 3 and 1 + 2 dimensions a constant magnetic field can 
spontaneously break chiral symmetry no matter how small the coupling constant G is, while in mag netic field free 
environment the chiral symmetry can be preserved when G is smaller than a critical value G c fl2l. fl3 . Il51 . 

From the magnetic catalysis effect, it is clear that the external magnetic field stimulates QCD condensation, with 
stronger magnetic held, the the value of QCD condensation or dynamic mass is bigger. On the other hand, as 
temperature participating in, the temperature weakens the condensation, with higher temperature, the value of QCD 
condensation or dynamic mass is smaller. The effect to condensation from magnetic held and temperature are opposite, 
hence putting these two conditions together can lead us to study their corporate contributions to the condensation. 
Beside that, in the chiral limit, when the temperature is high enough (exceeding a critical value T c ), the system shall 
undergo a phase transition, the broke chiral symmetry will restore, and now with the presentation of external magnetic 
held, T c should be relative to the magnetic held, hence one of our purpose to study the influence of magnetic held 
and temperature to QCD condensation is to identify the function relation of magnetic held and critical temperature. 

Susceptibilities and critical coefficients are important parameters to evaluate the chiral phase transition of QCD. In 
quantum held theory various susceptibilities are the linear responses of QCD condensation to various variables (e.g. 
temperature, current mass, chemical potential etc). In Ref. [ 161 ] . it is found that the nonlinear susceptibilities are 
correlated with the cumulant of baryon-number fluctuations in experiments, and in Ref. E3 ’s work, the researchers 
have calculated these susceptibilities without magnetic held and established the correlation, but in their work this 
correlation still has some deviation from the experiments, since in the relativistic heavy-ion collisions there are the 
presence of strong magnetic held, we expect the magnetic held could amend such deviation, hence in this paper we 
study the susceptibilities in NJL model with magnetic held in order to support further work. 
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The following of this paper is organized in such a way: In the section II of this paper, we give the numerical results 
of dynamical mass dependence on magnetic field and temperature. In section III we study the susceptibilities with 
non-zero current mass, compare them to the chiral limit case, and discuss their phase transition properties. In this 
section we also calculate the critical coefficients for different susceptibilities. In Appendix A, we have proposed a 
new method that is different from Schwinger proper time, its result is equivalent to other methods’, but when there 
is only external constant magnetic field rather than external electric field or external electromagnetic field, this new 
method is more convenient than the other method, at last of this appendix we discuss how this new method serves 
to thoroughly evaluate the contribution from infinitesimal imaginary term of the fermion propagator. In Appendix 
B, we give detailed deduction and discussion to the gap equation of two flavor NJL model with external magnetic 
field, because the electric charges of u d quarks are different, we need to prove the gap equation’s self-consistency 
with rigorous in this appendix. 


II. THE GAP EQUATIONS AND NUMERICAL RESULTS IN CHIRAL LIMIT 


In 1 + 3 dimensions, the bosonized two flavor NJL lagrangian with external magnetic field is 


C = ® Q — o — ij 5 ® 7? • t)i/j 



(1) 


n , q u 0 \ 2 1 

Q=\ Q l , Qu = g, Qd = qf=Qu,qd , 


( 2 ) 


£ 2 = a 2 + 7r 2 , tH = |7?| 2 = >V, 


2=1 


(3) 


T= (AW 3 ), T 1 = 


0 1 
1 0 


T 2 = 


0 -i 
i 0 


r 3 = 


1 0 
0 -1 


7T • T = 7 ToT 


(4) 


Here we have already assumed the current mass of fermion is zero, and A^ is the potential of external magnetic field, 
in order to serve the purpose of a constant external magnetic field, A M can be defined as 

{A 0 ,A 1 ,A 2 ,A 3 ) = (0, yz 2 ,-^ 1 ,!)), (5) 

x 1 and x 2 are the second and the the third component of time-space coordinates (x°, x 1 , x 2 , x 3 ). From Eq. the 
fermion propagator of different flavor with magnetic field should be 

= -, A f =idu. + q f eA fi , (6) 

- a 

throughout this paper, the pion condensation is thought to be zero (ir = 0), because its existence would violate parity, 
a detailed analysis is discussed in Appendix B. 

A more practical version of the fermion propagator is 


ijlf + a + a 

C P ) 2 - <J 2 ~ {Ilf) 2 ~ a 2 - qfeBa 12 ' 


(7) 


cr 12 =diag(l,-l,l,-l). (8) 

In Eq. ([7)1 . rather than employing Schwinger proper time method to rewrite the denominator of Sf, we propose a 
new method to deal with this denominator. The detailed deduction of this new method has shown in Appendix A, 
here we only give the final result. 

First of all, the general gap equation in two flavor NJL model with zero temperature is (the rigorous deduction of 
this equation is shown in Appendix B) 


cr 

G 



E Tr ^/- 

/ 


( 9 ) 
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Now we can employ the results from Appendix A, replace Eq. (IA19I) into Eq. ©, and simplify the gap equation, we 
have 




r+oo g —a s 


coth( | qf | eBs) ds , 


( 10 ) 


and by making a cutoff ^2 to the low limit of s’s integral, eventually we get the gap equation that is suitable for 
numerical calculation, 


47T 2 

~g = n ' 


f+oc g — a s 


d s + J2\qf\eB 
f 


coth(|g/|ei3s) — 


1 


\qf\eBs. 


ds. 


(ID 


In Eq. ED, Nf = 2, and the divergent part and convergent part are separated, hence the cutoff only affects the 
divergent part. 

At the finite temperature, the gap equation is slightly different with Eq. which is 


a 

G 


dr I d 4 x = -^Tr Sf, p = i 
J f 


( 12 ) 


and comparing with Eq. CD, the po operator has changed into Jy, with an appropriate eigenstate, there is 


— |m)o = *w m |m)o, = (2 to + 1)7tT, m e {■ ■ • ,-2,-1,0,1,2,• ■ ■}, 


then we can deduce an intermediate gap equation, 


4-7T 


+ OO +OO 


— = 2TY / \qf\eBldH 3 £ £ 
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w 2 j + 2n\qf\eB + II 2 + er 2 


(13) 


(14) 


m=— 00 n= 0 

In order to inherit the same cutoff ^ from Eq. (1A19I) . we choose to sum all of 2n\qf\eB up rather than add w m up, 

+00 +00 ^+00 


47T 2 

~G 


2 T^ 2 \qf\eB dU 3 ^ J ^(2 -S 0n ) e ~{u: m + 2 n\ qf \eB+n 3 +* )s ds 
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m =0 




■ coth(|g/|eBs) ds, 


(15) 


introducing the cutoff, 
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+00 / + 2 ? \ -a 2 s 

a \ 0 
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m —0 




coth( | qf | eBs) — 


\qf\eBs_ 


ds. 


(16) 


Now we can use Eq. ED to identify the relation between dynamical mass of NJL and constant external magnetic 
field eB , referring to Ref. j20|, for f n = 93MeV, = 138MeV and the current mass too = 5.5MeV, the value of 
cutoff A and G are 


A = 0.991GeV, G = 25.4GeV“ 2 . (17) 

Having set up these values, we are able to draw the a-eB relation in Fig. |TJ it is clear that when the strength of 
magnetic field increases, the dynamical mass will increases with it, the chiral symmetry is always broken, and when 
the magnetic field is strong enough, the dynamical mass has a nearly linear response to the magnetic field. 

Actually, when magnetic field is strong enough, the gap equation can be simplified to 

a _2 /*+oo as 

A V ds - (18) 

f A7 
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FIG. 1: The eB dependance of dynamical mass a when zero temperature. When eB = 0, cr = 0.25GeV. 
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FIG. 2: The eB dependance of critical temperature T c . When eB = 0, T c = 0.17GeV 


Because in Eq. m, when the cutoff 1/A 2 of integral variable s has made, there is 

lim coth(\qf\eBs) = 1, s > -F ( 19 ) 

eB—H- oo 

Of cause, Eq. (fl8l) can also be deduced from Eq. (IA16I) by throwing away the summation of n only leaving the n = 0 
term along. 

When temperature is not zero, it turns out there is a critical temperature T c , as long as the system’s temperature 
exceeds the critical temperature, the dynamical mass will be zero, which corresponds to the restoration of chiral 
symmetry. But for different magnetic field, the critical temperatures are different, as shown in Fig. [2] T c is increasing 
along with eB’s increase. While in Fig. [3] and Fig. [4] they give us magnetic field or temperature dependance of 
dynamical mass. In Fig. [3J the bigger the strength of eB 1 the upper the a-T curves, while in Fig. 01 the higher the 
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FIG. 3: The temperature dependance of dynamical mass a with different fixed eB. 
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FIG. 4: The eB dependance of dynamical mass a with different fixed temperatures. 


temperature, the lower the a-eB curves, and in this figure, the curve T = 0.2GeV is a little different with the other 
two, because in this situation, the magnetic field need to exceed a specific quantity to produce dynamical mass, while 
below this specific quantity, the chiral symmetry is unbroken. 

In Fig. [3] at different eB, the a-T curves are separating with each other, this a solid evidence that in NJL mean field 
approximation, there will have no inverse magnetic catalysis, and here we also can use a simple mathematics analysis 
to support this conclusion. Taking the cut off Eq. m for discussion, eB coth(\qf\eBs) is a monotone increasing 
function of eB, hence the RHS of Eq. (fl5l) is increasing along with eB increasing, this will obviously cause the a 
monotonously increasing whichever the temperature is. 


III. SUSCEPTIBILITIES AND CRITICAL COEFFICIENTS 

In this section, we refer to Ref. jH, HIJ’s works, use the known gap equation to study various susceptibilities. 
Because we do not consider the chemical potential in this article, there will be no susceptibilities related to chemical 
potential, such as quark number susceptibility Xq an d vector-scalar susceptibility Xvs- but after introducing external 
magnetic field in NJL model, we can define a new susceptibility related to eB, named as magnetic field susceptibility 
Xb- Therefore in this section, we are going to study three kinds of susceptibilities, xb-, chiral susceptibility Xc and 
thermal susceptibility xt- 

In order to study these susceptibilities, firstly, we need to enhance the NJL model Eq. m with current mass m, 
correspondingly, the new gap equation can be simply achieved by replacing a with a + m at the RHS of Eq. © and 
Eq. (1121) . We assume the cutoff and coupling constant is current mass irrelevant, hence the intermediate and final 
gap equation are 


47T 2 

~C!~ 


/*+oo / +oo 

a = 4v / 7rT(cr + m) \ eB / uC e 

f Jo V m=0 
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( 20 ) 
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/* + 00 / +°° 
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■ ds 


0 — (a+m) 2 s 
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coth(|g/|e5s) — 


\qf\eBs 


ds. 


( 21 ) 


Even with current mass, the properties of dynamic mass a are qualitatively similar to Fig. [l][4] 

Because {ipi/j) oc (— a), and we have assumed the coupling constant G and cutoff A are independent of current mass, 
temperature and magnetic field et cetera, therefore in this article, we define the susceptibilities as 


da da da 

Xc = dm’ XT = ~dT ’ XB = d{eB )' 


( 22 ) 


In Eq. (1211) . treating er as the implicit function of m, T and eB, a = a{m, T, eB), we can make partial differentiations 
of m, T and eB and get the corresponding equations for the susceptibilities. 
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FIG. 5: The temperature dependance of y c with different eB. 


For ensuring the equations are easy to read and analyse, firstly we define a few functions to represent some complicate 
formulas, 


/ +o° / +°Q \ 

Y e ~ U m s j 

TT n / 


g— (a+m) 2 s 


ds 


~£Z x m =0 

/»+oo / +°° 

+8 V^TY \Qf\eB 

f v m=0 


\ p-( CT + m ) 2 


Y e ~^ mS )e~ {a+rnrs V^ 


coth(|q/|ei?s) — 


\qf\eBs_ 


ds , (23) 


f T {a,m,T,eB) = 


/• + 00 / +°° 


8? T%N f T' z I ( ^(2m+l) 

•'a* x m=0 

/» + oo / + 00 

+87r*T 2 £> / |eB 

f v m=0 


„ a i e -( CT + m ) « 

2 ~-“'-' ' 6 ds 


Y (2m + l) 2 e"“- s ) e -(<r+m) 2 ^ 


coth(|g/|ei?s) — 


k/|eBs. 


ds, 

(24) 


/•+oo / +00 

fs(cr,m,T,eB) = 4 ^kTY Id/I / ( Y ' 

f 'm=0 


„ — (cr+m) 2 s 




[coth(|g/|eBs) — (|(?/|ei?s) csch 2 (|g/|e.Bs)] ds, (25) 


although these formulas are complicate, but they are all positive. 

Now by employing Eq. (ED, ©, HMD and (ED, we are able to deduce the equations for susceptibilities, 


1 


Xc = 


+ ^r(cr + m) 2 /rr 


- 1 , 


(26) 


er+m ' 47T 


XT = (1 + Xc) 


G , ^ , a 


G 

0r XT = - 4^2 (1+Xc) 


<9(RHS) 


<9T 


(27) 


Xb = ^2( 1 + Xc)(o’ + m)/ B . (28) 

For convenience of numerical calculation, Eq. (1271) provides a second formula to calculate thermal susceptibility, in this 
formula, (RHS) represents the RHS of Eq. (1211) . and the partial differentiation of (RHS) only operates on parameter 
T but treats a as a constant. 

The current mass is fixed to m = 5.5MeV, and the corresponding temperature or magnetic field dependance of 
susceptibilities are shown in Figs. ISlfTOl in Figs. [G] and [51 there seems to have anomalies, because the T = 0.2GeV 
curves in both figures have distinct values with T = O.OlGeV and T = 0.4GeV curves, this phenomenon is explainable, 
take Fig. [6] for example, making a cross-over analysis with Fig. [5J in Fig. [5] the peak of Xc is near T = 0.2GeV, 
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FIG. 6: The eB dependance of \c with different temperature. 
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FIG. 7: The temperature dependance of \t with different eB. 
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FIG. 8: The eB dependance of with different temperature. 



FIG. 9: The temperature dependance of xb with different eB. 
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FIG. 10: The eB dependance of \b with different temperature. 



T(GeV] 


FIG. 11: The temperature dependance of y c with different eB when m = 0. 


while at the both ends of T (T —> 0 and T —> 0.4GeV), \c approaches small values, therefore in Fig. [6] there comes 
the distinct difference between T = 0.2GeV and T = O.OlGeV T = 0.4GeV curves. In Figs. [Gl El EB the peaks of 
T = 0.2GeV curves tell us which value of eB is the crossover points. And in these figures, we can use them to identify 
the crossover properties, for example, in Fig. [Gl the T = 0.2GeV curve has no steep slopes around eB = 0.2GeV 2 , 
which means we can find out other crossover points at the nearby of T = 0.2GeV with magnetic field not far from 
eB = 0.2GeV 2 . 

From Figs. [5] [T] and 0 one can find the peaks are smooth, hence when the temperature increases from zero, the 
quark condensation will experience a crossover, and the crossover points is depending on external magnetic field, in 
these figures we can see when current mass has been considered, the crossover happens around T = 0.2GeV, besides, 
in Fig. [5] larger eB will make the T-\ c curve (take the peak as referent point) shift to right, in Fig. [71 larger eB 
not only makes T-\t curve shift to right but also makes the peak value increases, in Fig. [9l larger eB increases the 
values of but when eB is small, the crossover behavior of T-\b is not so obvious, therefore when magnetic field 
is weak, xb is not a good parameter to study crossover or phase transition behaviors. In these three figures, we can 
see that with the same magnetic field, all these crossover points in different susceptibilities are basically consistent, 
the consistency comes from the property of (1 + Xc)> because from Eqs. (1271) and (1751) . they are all (1 + Xc) -relevant. 

For comparison, we let current mass m be zero, and draw the corresponding \-T curves with different magnetic 
field, they are shown in Figs. [Til Fig. [T7] and Fig. [13] 

It is clear that when temperature approaches 0.2GeV, all of the three susceptibilities approaches infinity, from Eq. 
(E3l) and Eq. (ESI), one can find out a reasonable explanation, because when T crosses a critical point in the chiral 
limit, the dynamic mass a is zero, chiral symmetry is restored, this causes / m (0, 0, T, eB) being zero, which leads to 
the infinity of Xc■ For the other two susceptibilities xt and xb , although /t(0, a , T, eB) and /s(0, a, T, eB) approach 
zero, but apparently the approaching rate is lower than f m (0,a,T,eB), this ensures them also have infinity limits. 
After all, when m = 0, there is no necessary to study the susceptibilities’ properties at the T > T c area, because 
they are all infinities. These properties ensure us at chiral limit it is a second order phase transition. But one thing 
deserves mention here, at chiral limit, T = 0.2GeV is not the upper limit of phase transition, although from Figs. [11] 
Fig. IT71 and Fig. [13] it seems so, but from Fig. [2] we can see that bigger eB will linearly increase T c , eventually T c 
can exceed 0.2GeV with no upper limit. 

In order to study the second order phase transition, we also have studied the critical coefficients under chiral limit 
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T[GeV] 

FIG. 12: The temperature dependance of \t with different eB when m — 0. 
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FIG. 13: The temperature dependance of \b with different eB when m = 0. 


(m = 0), near the phase transition point (T is near and small than T c ), the definition of these coefficients are 

t=l -%r, (29) 

-L c 

a~t a , XT~t~ 1T , XB~t-ye. (30) 

a, 7 c, 7 t and 7 b are the critical coefficients that we are going to find out, for calculating a, on can draw the ln(<r)-ln(f) 
diagram and make a linear fitting, the slope of linear fitting line is a. The linear fittings are shown in Figs. 11411151 
[T6l fin they are the logarithmic relations between these critical coefficients and t, in these figures we have taken 
two different magnetic fields for comparison. In Table (T) we give the four kind of critical coefficients with different 
magnetic field. 



FIG. 14: The critical coefficient ln(cr)- ln(t) relation when m = 0, the slope is a. 
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FIG. 15: The critical coefficient ln(xc)-ln(t) relation when m = 0, the slope is 7 C . 



FIG. 16: The critical coefficient ln(xT)- ln(f) relation when m = 0, the slope is 7 t 



FIG. 17: The critical coefficient ln(\B)-ln(t) relation when m — 0, the slope is 7 b 


TABLE I: Critical coefficients with different eB 


eB 

a 

7= 

7 T 

7 B 

0.01 

0.47 

1.00 

0.58 

0.45 

0.05 

0.47 

1.00 

0.59 

0.44 

0.11 

0.45 

1.01 

0.63 

0.42 

0.21 

0.47 

1.02 

0.60 

0.45 

0.31 

0.47 

1.04 

0.60 

0.46 
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IV. SUMMARY AND CONCLUSIONS 

In this article, we have developed a new method to study the NJL gap equation with external magnetic field, this 
method is consistent with other method, and it is also compatible with more complicate equations such as Dyson- 
Schwinger equation, beside, this method is particularly convenient for dealing with constant magnetic field problems, 
while for a general external electromagnetic field the proper time method is still a better way. In our work, we have 
successfully repeated other researchers’ works, these results are shown in Fig. [THU from these figures we know the 
’Magnetic Catalysis’ is not only enhance the QCD condensation but also increase the value of critical temperature. 
This is reasonable, because phenomenologically the external magnetic field restricts the movement modes of charge 
particles, gives more chance for any two particles to pair up, hence the QCD condensation increases. We have also 
calculated various susceptibilities, and all these susceptibilities imply that in chiral limit, the condensation has the 
second order chiral phase transition at finite temperature, while in non chiral limit, even the current mass is tiny, the 
condensation only experiences crossover. In this article we do not study chemical potential, because with chemical 
potential, there is a problem of chemical potential feed back, more than that, the presence of magnetic field will 
entangle the feed back, one can refer to the appendix for a simple impression, rigorous study of this problem is left in 
proceeding articles. In the chiral limit, the Wigner solution is trivial (a = 0), the susceptibilities of Wigner solution 
is also trivial \c = 0 (in Eq. (EH1) . let a —> 0 first, ant then let to -A 0, if one changes the limitation sequence, the 
result is infinity), hence we do not discuss in this paper. 


Appendix A: A New Method Different from Schwinger Proper Time 

For simplicity, Il£ is wrote as II M in this appendix. Unlike free fermion propagator, II i and II 2 are not commutable, 
hence it is not possible to find a representation in which all four II^’s eigenstates exist simultaneously, at least the 
eigenstates of IIi and II 2 can not exist simultaneously. But for the denominator of fermion propagator, we do not 
need to find the eigenstates for all 11^, instead we turn to find the eigenstates of II 2 . 

Temporarily, we define the eigenstate of II 2 as |n^,IIo,n 3 ), in which 11^ represents the eigenvalue of (Ilf + Ilf), 
therefore II 2 |II^, IIo, II 3 ) = (Ilf — — Ilf )|IIj_, Ifo, II 3 ). From the definition, we can see that the eigenstate only have 
three degrees of freedom, but to completely describe a particle there needs four degrees of freedom, therefore, in the 
following article, we are going to find out the rest degree of freedom. Besides, the operators IIo and II 3 are actually 
Po and P 3 , so the eigenvalues IIo and II 3 are equivalent to po and p$. 

For the eigenvalue of IIo and II 3 , their values vary continuously from —00 to + 00 , while the property of IIj_ 
are not identified yet. Because the operators IIj_, IIo and II 3 are commutable with each other, we can rewrite the 
eigenstate in a form of tensor product, |Ilf_, n 0 , II 3 ) = |n 0 )o 0 |n^)i2 0 |n 3 ) 3 , which allows us to study the property 
of U 2 ± separately. In the coordinates representation, |n^)i 2 lies in the lx 1 )! 0 |x 2 }2 space, a general idea is to find the 
expression of 12 (x 1 , x 2 |Ilf_)i 2 through differential equation, but here we are not going to directly deduce the expression 
for l 2 (x 1 ,X 2 |Ilf_)i 2 , instead we try to find an intermediate expression by introducing a new representation beyond the 
coordinates representation and momentum representation. 

In the two dimension Hilbert space which in the language of coordinate representation is the tensor space on the basis 
of lx 1 ): 0 |x 2 ) 2 , we can reorganize the basis tensor through eigenstates of ITi. Because II 1 = Pi - eAi = —p 1 + ^l^-x 2 , 
and the momentum operator p 1 and coordinate operator x 2 operate on different Hilbert spaces (lx 1 )! and |x 2 }2 
separately), then we can define (ITi),j ? 1 as the complete operators that cover the whole lx 1 )! 0 |x 2 )2 space, their 
eigenstates is 


l n LP) = y Iq^eB ^ 1 ® l a ’ ni ) 2 ’ XUl = + ^ ( Al ) 

Here the eigenstate |p)i is for p 1 rather than pi, p 1 ^)! = p\p)i, Pi\p)i = —p\p)i, one must caution the sign problem, 
actually in the context of this article, all the eigentstates of coordinates and momentum are eignestates of x M or pfo 

In Eq. m, the factor \ q ^ eB is used to normalize the basis tensors, hence we have the normalized, complete and 
orthogonal relations, 


(n , 1 ,p'\iL 1 ,p} = 6(n-n')6(p-p'), 


(A2) 


dnidp|ni,p)(ni,p| = 101, 


(A3) 
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and by employing the equations p 2 \x 2 )2 = \x 2 )2 and ^ 1 |p 1 )i = — ’i-^r\p 1 )i, we have 


o \ Q 

n 2 |ni, P ) = -iq f eB[ — - )|ni,p). 


(A4) 


Now with this preparation of eigenstates |IIi,p), we are able to establish a differential equation for |IIj_) 12 - For 
operator fl^. there is 


<n 1 ,p|ni|i^.) = ni<n 1 ) p|i^) > 

while ft]/ = ftf + ft|, with Eq. (1A4I) . there is also 


(A5) 


(n^piniini) = n^n^ni) - q 2 e 2 B 2 




2 

<ni,p|ni>, 


(A6) 


put Eq. (IA5I) and Eq. (IA6I) together, we have the differential equation of (IIi,p|n[/), the variables are III and p , 


ni<n„ P |ni> = (n; - )(n„p|ni). 


If we assume (IIi,p|n5_) has the following form, 


(ni,p|ni) = * = / ^n„ 

a is an arbitrary complex constant, then Eq. (1X71) can be simplified to 

d 2 h(z) f 1 z 2 \ , n/ 1 

— + { n+ 2-T) h{z} = l> ’ n= 2WlB*2' 


(A7) 


(A8) 


(A9) 


Eq. (IA9D corresponds to Weber differential equation, the solutions of this differential equation depend on the property 
of n, here for the physical purpose, we need the solutions be convergent when z —> ±oo, therefore this requirement 
constrains n be a nonnegative integer, which also determines the property of 11/ as 11/ = (2 n + l)\qf\eB, n € 
{0,1,2,---}. In the following article, we rewrite function h(z) as h n (z) to show the n-dependance of the solution. 
The general expression of h n (z) is 


h n (z) = e 4 z r 


[n/2] , n . / l—n \ 

l 2> k ^ ~> k 

k -0 


k\ 



(A10) 


where [n/2] is a nearest integer function, and (—§)fc, k obey the following rule, for an arbitrary real number A, 

a positive integer fc, there is 


(A)o — 1, (A )k — — (A + k — 1)(A + k — 2) • • • (A + 1)A. (All) 

Now we have identified II]/ and h n (^J i), let’s go back to Eq. (IA 8 I) to identify parameters a and c. For the 

same physical consideration, we need (Il^pln/) be convergent when p, Dp —> ±oo, this leaves a no other choice but 
a pure imaginary number, so we redefine a as ia , now a is an arbitrary real number, it represents the a hidden degree 
of freedom, hence in the following discussion we correct the eigenstate |n/)i 2 to \n, a). 

c is a normalization factor, to identify c’s value, we shall give the orthogonal and complete relations of h n (z ) first, 


J h m (z)h n (z) dz = n\V2Tr5 mn , 


(A12) 


S(x 


+oo J 

y) = V] , MaOMy)- 


n —0 


(A13) 
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In order to normalize \n,a), it is reasonable to assume c is n-dependant, let c — > c n , 

1 


(IIi,p|n, a) =c n e la(p+ 2l) h n ( x 


jIIi)) c n — 


k/|ei? ’ \n\2i: ^/\qJ\eBn 

Therefore eigenstate | n,a) has the following normalized and orthogonal relation, 

(to, a'\n, a) = J dUidp (to, a / |ni,p)(IIi,p|n, a) = 6 mn S(a — a'). 


(A14) 


(A15) 


A short summary, for operators (ft^, li 0 , II 3 ), we can use their eigenstates |n 0 ,II 3 ;n,a} = |n 0 )o <8> |ro, < 2)12 <8> 1113)3 
as a set of complete basis tensors in four dimension Hilbert space. 

Normally, we need to calculate Tr.S', for example in Eq. <[9|), the operator Tr is representation irrelevant, hence with 
constant external magnetic field, Tr S can be treated as 


/f +OO 

dUodlls / da^^( lip,II 3 ; n, a\ tr 5 |IIq, II 3 ; n, a) 

n—0 


+00 


— 2 a dn 0 dn 3 (n 0 , n 3 |no, n 3 ) 


(2 - Son) 


n—0 


— 2n\qf\eB — n| — cr 2 


da (n, a\n, a). 


(A16) 


In order to cancel J d 4 x in LHS of Eq. ©, J d i x must be extracted from Eq. (IA16I) . we can employ the following 
relations to serve this purpose, 


(n 0 ,n 3 n 0 ,n 3 ) — 2 j 

^ dxodxo, 

(A17) 

J da (n, a\n, a) - |g/|e ^ 

j " dx\dx2, 

(A18) 


replacing them into Eq. (IA16I) . and let n 0 have a Wick rotation, it will be simplified to 


.\qf\eBa 

[ d A x 

j dU 0 dU 3 t-~v (2 — 5on) 


% n J 

1 ( 2 H ) 2 ^ o H 2 0 + 2n\q f \eB + Ul + a^ 


,\q f \eBa 

[ d 4 x 

[ S 2 M / + °° ds 

' ' ’ n=0 -'° 


,\q f \eBa 
47T 2 J 

( d 4 x 

r +00 e -a 2 s 

/ -coth( \qf \eBs) ds. 

1 0 s 

(A19) 


In this appendix, we have developed a new method that is different from Schwinger proper time method. The 
original idea of developing this new method is to evaluate how much influence does the infinitesimal imaginary term 
of fermion propagator has to the final results. Because in quantum field theory, the infinitesimal imaginary term 
(which always presents in the form of ie) does not always play a role as pointer of integral path, in some case it can 
cause remarkable adjustment to the calculations, for example, in the case of zero temperature and finite chemical 
potential p, in Minkowski space (lol l30j|. the fermion propagator is 

s ^^=T2 — 2 ~ h+ ,r ^ —r’ (A2o) 

whose infinitesimal imaginary term is /r-dependent, then after the Wick rotation of po, there will be nonzero residue 
counted in. Out of the consideration that the influence of a external magnetic field might affect infinitesimal imaginary 
term, therefor a thorough method is needed. For a free fermion propagator without any external field, the infinitesimal 
imaginary term is ie, this is well-known in quantum field theory, but when external fields or external elements interfere, 
it is not safe to claim that the imaginary term is ie, a convincible example is the case that the system has finite 
chemical potential, in Minkowski space, the infinitesimal imaginary term is total chemical potential dependance, 
therefore one need to properly deduce the the imaginary term with caution, generally there have two methods to 
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deduce the fermion propagator with imaginary term, one is canonical quantization, the other is path integral. Here 
we introduced a convenient trick for path integral I18l Il9j | , the definition of partition function of quantum field theory 
is 


Z = 


Texp 




(A21) 


Now introducing a factor (1 — irj) to change the expression of partition function, which will derive a 77 -dependent 
Lagrangian, 


Z = lim ( 0 
?j—> 0+ 

= lim 

rj— >0+ 


^0 T exp | — i{l — irj) J H ctt 


dip (ip 


T exp | — i(l — 777 ) 

= lim [ D'lpD'tp e® f dx Cri . 
v^o+ J 


p + OO 


H dt. 


iP 


(A22) 


In NJL model with external magnetic field, comparing with Eq. the improved fermion propagator is 

P + a 


S = 


QJ ) 2 — a 2 + 70 ( 77 ) ’ 


0 ( 77 ) = 7?(|nf + er — eBa 42 ). 


(A23) 


From previous discussion, we can see that is quantized by magnetic field to (2 n+ 1 )eB, hence the imaginary term 
0 ( 77 ) is permanently positive, it is equivalent to ie , the influence is trivial. But for the finite chemical potential case, 
one can prove that the tricks do give the right fermion propagator with 77 . 

So far, we have used this new method to deduce the gap equations such in NJL model with magnetic field. In fact 
these gap equations are the same as the ones deduced through Schwinger proper time method, the difference is the 
process, as mentioned before, this new method is developed to evaluate the infinitesimal imaginary term of fermion 
propagator, and of cause, the influence from imaginary term is trivial, but anyway, it deserves a mention how this been 
proved. This method is very convenient to study NJL problems with external magnetic field, because one does not 
need to find expression for (x\S\y) at the first place, and it is also qualified for further study such as Schwinger-Dyson 
equations, but unfortunately this method can not simplify the complexity neither when higher order contributions 
such as fermion loop are considered. 


Appendix B: Deduction and Discussion of two flavor NJL gap equation with magnetic field 
With Eq. (Ill), we are able to write the free energy as 

T = ^(cr 2 + 7 T 2 ) + N c iTif In \ij) + ef ® Q - (a + 77 s <g) if • f)], 

Z(j 


(Bl) 


the lower index ’/’ of Tr/ means beside the trace of spinor space and integral of 4 dimension coordinate or momentum 
space there also has a trace of flavor space. 

From Eq. m one can deduce the gap equation for 7?, 


7T 

G 


d 4 x = i Tr/Kfy 5 Cg) t)S]. 


(B2) 


When = 0, one can easily prove that n can have trivial solution 7 ? = 0, normally in NJL model, this solution is 
widely accepted, and most of the researches are study the gap equation of a. But in our case here, there are two 
reasons that we can not directly employ this conclusion, firstly, from Eq. 0 one can find in the denominator of Sf, 
the non-commutable relation of n{ and n| produces the qfeBa 12 term in the spinor space, secondly, the difference 
of electric charges of u d quarks produces a non-identity matrix Q in the flavor space, therefore it is hard to deduce 
the final gap equations for a and if, and 7? might not have a trivial solution. In this appendix, we are going to prove 
that 7 ? still can have a trivial solution. 

The complete propagator S in Eq. 


is 


S = 


1 


i + £g) Q — (a + 77 s <8> 7? • r) 


(B3) 
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In this propagator, because of the existence of <8> Q, it is not easy to construct simple quadratic terms in the 
denominator. 

Since our purpose is to prove if can have trivial solution, we assume it is, then replace 7? = 0 into Eq. (|B2I) . if the 
LHS and RHS of Eq. (1B2I) are equal, our assumption is proved. 

When 7? = 0, the LHS of Eq. (IB2D is 0, now we are going to prove that the RHS of Eq. (IB2D is also 0. For the 
propagator, there is 


S | ff=0 = 


_ = [ 

i + <8> Q — <J l O 


O 


i/i d -o 


S u O 

o s d 


Therefor, we have 


(iq 5 <g> t 3 ){S |*=o) = 


n 5 s u O \ 
o -i 7 5 S d ) ’ 


replacing it into the RHS of Eq. 

Ti 7 [(7 7 5 ® r 3 )(S| ff=0 )] = Tr(n 5 £„) - Tr(i 7 5 S d ) = 0. 
As for r 1 ’ 2 of r, there definitely has 

Tr/[(i 7 5 ® r 1 ’ 2 )(S'| s=0 )] = 0. 


(B4) 


(B5) 


(B6) 


(B7) 


Now we have proved that with the presence of external magnetic field, we still can find trivial solution for if. In fact 
this is rational, because iy 5 ® 7?t in the denominator of complete fermion propagator violates parity, hence in many 
papers if is set as 0 [31[ 153 ■ 

Through this conclusion, now we are able to simplify the gap equation of <7, 


§ Jd 4x = iTr f (S) =iTr f (Sy= 0 ) = *£>(S/). 


(B8) 


Generally speaking, a should be a (4 x 4) ® (2 x 2) matrix, but from previous discussion we know that the (2 x 2) part 
(flavor space) is eventually simplified by the Try operator and 7 ? = 0, if one demands if 7 ^ 0, then the gap equations 
of a and if will definitely not be such simple. On the other hand the (4 x 4) part (spinor space) is simplified to scaler 
at the beginning in Eq. ©, because <r is introduced here as an auxiliary scaler field. 

Actually in NJL model, it is not obvious that the dynamic mass E should be a (4 x 4) matrix (here for convenience 
we only consider one flavor NJL model, the flavor space is neglected), because either mean field approximation or 
auxiliary field method, it provides only a scaler rather than a matrix to the dynamic mass (or non-approximate 


fermion self energy). While in Ref. [33]’s work, the authors had proposed a method to find out the matrix form of self 


energy, they used Fierz transformation to identify the matrix structure. But in a matter of fact, the result is trivial 
as long as chemical potential is not involved. 

Instead of employing the method used in Ref. J33j, here we use Dyson-Schwiger equation with contact interaction 
model to clarify this problem, the Dyson-Schwinger equation of fermion self energy is 


-iT,(x, y) = g 2 J dz 1 dz 2 YD^x, z±)S(x, z 2 )T ,/ (z 1 -, z 2 , y). (B9) 

For the contact interaction model, the gluon propagator and dressed vertex are simplified into 

D flu (x,z 1 ) = Xg flv S(x - zi), T v (zr,z 2 ,y) = 7 1/ <5(zi - y)S(z 2 - £ 1 ), (B10) 

replacing them into Eq. (IB9I) . 

E = i\g 2 (x\^ * 7m 1j;). (Bll) 

5 0 E 

In Eq. (1B11|) . it is rather than TrS 1 in the NJL gap equation, this difference naturally leads to the assumption 

that E is a 4 x 4 spinor space matrix as long as it complies with some specific symmetry. Normally E should be 
E = a + there is no 7 5 and 7 5 7 ^ terms because they violate parity, there is also no a^ terms because if there 
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was, in the LHS of Eq. (IB11I) . the a^ terms of E can not find corresponding terms in the RHS because 7 / V T 1 / 7 At = 0. 
On the other hand, a „ in E can be absorbed by hence their contributions are trivial, this leaves only a contributes 
to the self energy E, therefore a E can be simplified to a scaler. 

In this article, although magnetic field has been involved, we can also prove that a^ in E = a M 7^ + <7 is trivial, and 
a scaler E is the simplest and self-consistence solution for our magnetic field case. Assuming E = a M 7 M + a (there 
still not has 7 s , 7 5 7 M , a terms for the same reason as mentioned above), from Appendix A, we know 

[lx l-f-lx) d A x = a /4 + ba 12 . (B12) 

J P - «m7 m - 0 

And for the spinor matrices, there is 

7 ^ 7 ,, = 0, (B13) 

hence there is 

,/^ = U 9 v(/<* |S W ^) 7 „ = 4U 9 V (B14, 

the LHS and the RHS are both scalers, the assumption of scaler E is reasonable. 
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